Proceedings of the National Academy of Sciences

Plant cuticles are composed of wax and cutin, and evolved in the land plants as a hydrophobic boundary that reduces water loss from the plant epidermis. The expanding maize adult leaf displays a dynamic, proximodistal gradient of cuticle development, from the leaf base to the tip. Laser microdissection RNA Sequencing (LM-RNAseq) was performed along this proximodistal gradient, and complementary network analyses identified potential regulators of cuticle biosynthesis and deposition. Correlations between cuticle development and cell wall biosynthesis processes were identified, as well as evidence of roles for auxin and brassinosteroids. In addition, our network analyses suggested a previously undescribed function for PHYTOCHROME-mediated light signaling during cuticular wax deposition. Genetic analyses reveal that the phyB1 phyB2 double mutant of maize exhibits abnormal cuticle composition, supporting predictions of our coexpression analyses. Reverse genetic analyses also show that phy mutants of the moss Physcomitrella patens exhibit abnormal cuticle composition, suggesting a role for light-stimulated development of cuticular waxes during plant evolution.

Understanding how climate refugia and migration over great distances have facilitated species survival during periods of past climate change is crucial for evaluating contemporary threats to biodiversity. In addition to tracking a changing climate, extant species must face complex, anthropogenically fragmented landscapes. The dominant conifer species in the mesic temperate forests of the Pacific Northwest are split by the arid rain-shadow of the Cascade Range into coastal and interior distributions, with continued debate over the origins of the interior populations. If the Last Glacial Maximum extirpated populations in the interior then postglacial migration across the arid divide would have been necessary to create the current distribution, whereas interior refugial persistence could have locally repopulated the disjunction. These alternative scenarios have significant implications for the postglacial development of the Pacific Northwest mesic forests and the impact of dispersal barriers during periods of climate change. Here we use genotyping-by-sequencing (ddRADseq) and phylogeographical modeling to show that the postglacial expansion of both mountain hemlock and western redcedar consisted largely of long-distance spread inland in the direction of dominant winds, with limited expansion from an interior redcedar refugium. Our results for these two key mesic conifers, along with fossil pollen data, address the longstanding question on the development of the Pacific Northwest mesic forests and contrast with many recent studies emphasizing the role of cryptic refugia in colonizing modern species ranges. Statement of SignificanceUnderstanding whether habitat fragmentation hinders range shifts as species track a changing climate presents a pressing challenge for biologists. Species with disjunct distributions provide a natural laboratory for studying the effects of fragmentation during past periods of climate change. We find that dispersal across a 50-200-km inhospitable barrier characterized the expansion of two conifer species since the last ice age. The importance of migration, and minimal contribution of more local glacial refugia, contrasts with many recent studies emphasizing the role of microrefugia in populating modern species distributions. Our results address a longstanding question on the development of the disjunct mesic conifer forests of the Pacific Northwest and offer new insights into the spatiotemporal patterns of refugial populations and postglacial vegetation development previously unresolved despite decades of paleoecological studies.

By masterfully balancing directed growth and passive mechanics, plant roots are remarkably capable of navigating complex heterogeneous environments to find resources. Here we present a theoretical and numerical framework which allows us to interrogate and simulate the mechanical impact of solid interfaces on the growth pattern of plant organs. We focus on the well-known waving, coiling and skewing patterns exhibited by roots of Arabidopsis thaliana when grown on inclined surfaces, serving as a minimal model of the intricate interplay with solid substrates. By modelling growing slender organs as Cosserat rods that mechanically interact with the environment, our simulations verify hypotheses of waving and coiling arising from the combination of active gravitropism and passive root-plane responses. Skewing is instead related to intrinsic twist due to cell file rotation. Numerical investigations are outfitted with an analytical framework that consistently relates transitions between straight, waving, coiling and skewing patterns with substrate tilt angle. Simulations are found to corroborate theory and recapitulate a host of reported experimental observations, thus providing a systematic approach for studying in silico plant organs behavior in relation to their environment. SignificancePlant roots exhibit an exceptional ability to navigate in heterogeneous soil environments while overcoming obstacles. Our study combines theory and experimental observations to interrogate and simulate the mechanical impact of obstacles on organ growth. As a test case we focus on well-known observations of waving, coiling and skewing growth patterns of Arabidopsis thaliana roots grown on inclined substrates. Overall, our study explains a broad set of experimental observations through the minimal ingredients of gravitropism and passive mechanics. Our numerical framework provides an in silico laboratory, yielding quantitative insight into the dynamics of growing organs at the intersection of active processes and passive mechanics, applicable beyond plants to any slender growing system, from neurons or fungal hyphae to novel soft robots.

In plants, host-pathogen coevolution often manifests in reciprocal, adaptive genetic changes through variations in host nucleotide-binding leucine-rich repeat immune receptors (NLR) and virulence-promoting pathogen effectors. In grass powdery mildew (PM) fungi, an extreme expansion of a RNase-like effector family, termed RALPH, dominates the effector repertoire, with some members recognized as avirulence (AVR) effectors by cereal NLR receptors. We report the structures of the sequence-unrelated barley PM effectors AVRA6, AVRA7 and allelic AVRA10/AVRA22 variants, which are detected by highly sequence-related barley NLRs MLA6, MLA7, MLA10, and MLA22, and of wheat PM AVRPM2 detected by the unrelated wheat NLR PM2. The AVR effectors adopt a common scaffold, which is shared with the ribonuclease (RNase) T1/F1-family. We found striking variations in the number, position, and length of individual structural elements between RALPH AVRs, which is associated with a differentiation of RALPH effector subfamilies. We show that all RALPH AVRs tested have lost nuclease and synthetase activities of the RNase T1/F1- family and lack significant binding to RNA, implying that their virulence activities are associated with neo-functionalization events. Structure-guided mutagenesis identified six AVRA6 residues that are sufficient to turn a sequence-diverged member of the same RALPH subfamily into an effector specifically detected by MLA6. Similar structure-guided information for AVRA10 and AVRA22 indicates that MLA receptors detect largely distinct effector surface patches. Thus, coupling of sequence and structural polymorphisms within the RALPH scaffold of PMs facilitated escape from NLR recognition and potential acquisition of diverse virulence functions.

Thylakoid-free cyanobacteria are thought to preserve ancestral traits of early-evolving organisms capable of oxygenic photosynthesis. However, and until recently, photosynthesis studies in thylakoid-free cyanobacteria were only possible in the model strain Gloeobacter violaceus. Here, we report the isolation, biochemical characterization, cryo-EM structure, and phylogenetic analysis of photosystem I from a newly-discovered thylakoid-free cyanobacterium, Anthocerotibacter panamensis, a distant relative of the genus Gloeobacter. We find that A. panamensis photosystem I exhibits a distinct carotenoid composition and has one conserved low-energy chlorophyll site, which was lost in G. violaceus. These features explain the capacity of A. panamensis to grow under high light intensity, unlike other Gloeobacteria. Furthermore, we find that, while at the sequence level photosystem I in thylakoid-free cyanobacteria has changed to a degree comparable to that of other strains, its subunit composition and oligomeric form might be identical to that of the most recent common ancestor of cyanobacteria.

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) is foundational to life on Earth, catalyzing carbon dioxide (CO2) fixation to generate biomass. However, Rubisco is a slow and inefficient enzyme that has proven challenging to engineer. We applied the structure-informed machine learning (ML) model ESM-IF1 to identify plausible amino acid sites in the large subunit of Nicotiana tabacum Rubisco to target for directed evolution. ML-assisted library design followed by selection in Rubisco-dependent Escherichia coli identified multiple enriched variants displaying improved catalytic efficiency. Several improved variants carried amino acid changes not found in the evolutionary lineage of plants, despite being assembly competent in plant chloroplasts, demonstrating that ML-assisted protein design can explore functional sequence space beyond what is observed from natural sequence diversity. Most prominently, the T391I substitution improved carboxylation rate by 29% and aerobic carboxylation efficiency by 43%. Our findings demonstrate the utility of ML-assisted evolution for engineering Rubisco with improved carboxylation efficiency and potential for enhancing crop productivity.

Aquatic life is strongly structured by light gradients, with gradual decrease in light intensity and differential attenuation of sunlight wavelengths with depth. How phytoplankton perceive these variations is unknown. By providing the first in vivo quantitative assessment of the action of marine diatom phytochrome photoreceptors (DPH), we show that they efficiently trigger photoreversible responses across the entire light spectrum, unlike current models of phytochrome photosensing. The distribution and activity of DPHs in the environment indicate that they are extremely sensitive detectors of spectral light variations related to depth and optical properties of the water column in temperate and polar oceans, revealing a completely novel view of how light is perceived in the marine environment.

Mounting evidence suggests that plants engage complex computational processes to quantify and integrate sensory information over time, enabling remarkable adaptive growth strategies. However, quantitative understanding of these computational processes is limited. We report experiments probing the dependence of gravitropic responses of wheat coleoptiles on previous stimuli. First, building on a mathematical model that identifies this dependence as a form of memory, or a filter, we use experimental observations to reveal the mathematical principles of how coleoptiles integrate multiple stimuli over time. Next, we perform two-stimulus experiments, informed by model predictions, to reveal fundamental computational processes. We quantitatively show that coleoptiles respond not only to sums but also to differences between stimuli over different timescales, constituting first evidence that plants can compare stimuli - crucial for search and regulation processes. These timescales also coincide with oscillations observed in gravitropic responses of wheat coleoptiles, suggesting shoots may combine memory and movement in order to enhance posture control and sensing capabilities.

Witchweeds (Striga spp.) are parasitic plants that severely constrain cereal production across sub-Saharan Africa, threatening food security for millions of people (Runo and Kuria, 2018). Striga infection begins when dormant seeds germinate in response to host-derived biomolecules, primarily strigolactones which plants emit to regulate shoot branching and to communicate with beneficial microbes. This obligate dependence on host signals can be exploited for Striga control through suicidal germination, whereby strigolactone-like compounds induce parasite germination in the absence of a host. Although this strategy proved highly effective during Striga eradication efforts in the United States using ethylene gas as a Striga germination inducer (Eplee, 1975; Iverson et al., 2011), its deployment in Africa has been limited by capacity to synthesize cost-effective strigolactone-like Striga germination inducers. Here, we show that structure-guided in silico screening of chemical libraries using AlphaFold2-modeled receptor-ligand interactions improve the efficiency and likelihood of identifying previously unknown strigolactone analogs. Using this approach, we identify a structurally simple synthetic lactone scaffold that induces Striga germination at nanomolar concentrations. These results present new avenues for the development of strigolactone analogs and support revisiting suicidal germination as a practical Striga control strategy in Africa.

The size and density of stomatal pores limit the maximum rate of leaf carbon gain and water loss (gmax) in land plants. The limits of gmax due to anatomy, and its constraint by the negative correlation of stomatal size and density at broad phylogenetic scales, has been unclear and controversial. The prevailing hypothesis posits that adaptation to higher gmax is typically constrained by geometry and/or an economic need to reduce the allocation of epidermal area to stomata (stomatal-area minimization), and this would require the evolution of greater numbers of smaller stomata. Another view, supported by the data, is that across plant diversity, epidermal area allocated to guard cells versus other cells can be optimized without major trade-offs, and higher gmax would typically be achieved with a higher allocation of epidermal area to stomata (stomatal-area increase). We tested these hypotheses by comparing their predictions for the structure of the covariance of stomatal size and density across species, applying macroevolutionary models and phylogenetic regression to data for 2408 species of angiosperms, gymnosperms, and ferns from forests worldwide. The observed stomatal size-density scaling and covariance supported the stomatal-area increase hypothesis for high gmax. A higher gmax involves construction costs and maintenance costs that should be considered in models assessing optimal stomatal conductance for predictions of water use, photosynthesis, and water-use efficiency as influences on crop productivity or in Earth System models.

Reactive oxygen species (ROS) have been implicated repeatedly in multiple signaling processes in plants but the underlying mechanisms and roles remain enigmatic. Here, we developed live imaging of apoplastic ROS at the root surface. Different signals, including auxin, extracellular ATP and RALF1 peptide, all induce cytosolic calcium transients and apoplastic ROS bursts. Genetic and optogenetic manipulations identified calcium transients as necessary and sufficient for ROS bursts via activation of NADPH oxidases RBOHC and RBOHF. Apoplastic ROS bursts are not required but rather limit the gravity-induced root bending. Root bending is sensed by stretch-activated calcium channel MCA1 leading to NADPH oxidase activation at the stretched side. The resulting ROS production stiffens cell wall for better soil penetration. Apoplastic ROS thus provides a means to balance tissue flexibility and stiffness to efficiently navigate soil.

The evolutionary history of plant interactions with necrotrophic pathogens that feed on dying host cells and their virulence mechanisms remains fragmentary. We have isolated the barley gene Scs6, which is required for the necrotrophic fungus Bipolaris sorokiniana isolate ND90Pr to cause spot blotch disease. Scs6 is located at the disease resistance gene locus Mildew locus a (Mla) and encodes an intracellular nucleotide-binding leucine-rich repeat receptor (NLR). In transgenic barley, Scs6 is sufficient to confer susceptibility to ND90Pr in accessions naturally lacking the receptor, resulting in infection-associated host cell death. Expression of Scs6 in evolutionarily distant Nicotiana benthamiana reconstitutes a cell death response to an uncharacterized non-ribosomal peptide effector produced by ND90Pr-specific non-ribosomal peptide synthetases (NRPSs) encoded at the VHv1 virulence locus. Our data suggest that the heat-resistant effector directly activates the SCS6 receptor. Scs6 is an allelic variant of functionally diversified Mla resistance genes each conferring strain-specific immunity to barley powdery mildew isolates with a matching proteinaceous pathogen effector. Domain swaps between MLA and SCS6 NLRs and expression of the resulting hybrid proteins in N. benthamiana reveal that the SCS6 leucine-rich repeat domain is a specificity determinant for the NRPS-derived effector to activate the receptor. Scs6 evolved after the divergence of barley from wheat and is maintained in several wild barley populations with an incidence of 8%, suggesting a beneficial function for the host. Evolution of the bona fide immune receptor SCS6 targeted by the NRPS-derived effector was key for the emergence of strain-specific spot blotch disease in domesticated barley.

Plant height is an important agronomic trait with a significant impact on grain yield, as demonstrated by the positive effect of the REDUCED HEIGHT (RHT) dwarfing alleles (Rht1b) on lodging and harvest index in the "Green Revolution" wheat varieties. However, these gibberellic acid (GA) insensitive alleles also reduce coleoptile length, biomass production, and yield potential in some environments, triggering the search for alternative GA-sensitive dwarfing genes. Here we report the identification, validation and characterization of the gene underlying the GA-sensitive dwarfing locus RHT25 in wheat. This gene, designated as PLATZ-A1 (TraesCS6A02G156600), is expressed mainly in the elongating stem and developing spike and encodes a plant-specific AT-rich sequence- and zinc-binding protein (PLATZ). Natural and induced loss-of-function mutations in PLATZ-A1 reduce plant height and its over-expression increases it, demonstrating that PLATZ-A1 is the causative gene of RHT25. PLATZ-A1 interacts physically and genetically with RHT1 (DELLA), and both genes have stronger effects on plant height in the presence of the wildtype than in the presence of the mutant allele of the other gene. These results suggest that PLATZ1 can modulate the effect of DELLA on wheat plant height. We identified four natural truncation mutations and one promoter insertion in PLATZ-A1 that are more frequent in modern varieties than in landraces, suggesting positive selection during wheat breeding. These mutations can be used to fine-tune wheat plant height and, in combination with other GA-sensitive dwarfing genes, to replace the GA-insensitive Rht1b alleles to search for grain yield improvements beyond those of the Green Revolution varieties. Significance StatementWe have identified and characterized a previously unknown gene controlling plant height in wheat and named it PLATZ1. Mutations in PLATZ1 reduce plant height while its overexpression results in taller plants. PLATZ1 is expressed mainly in elongating stems and developing spikes and interacts physically and genetically with the "Green Revolution" dwarfing gene REDUCED HEIGHT 1 (RHT1). We discovered five natural mutants in the A genome copy of PLATZ1 in common wheat that have been favored during breeding, suggesting an overall positive effect on wheat performance. These mutations can be used to fine-tune wheat plant height and, eventually, to replace the RHT1 dwarfing alleles that impose limitations on planting depth and grain yield potential in some environments.

Chlamydomonas acclimates to repeated low (LL) or high light (HL) days by changing the abundance of photosynthetic complexes and the ultrastructure of its thylakoid membranes. These phenotypes persist through the night phases, suggesting a readiness for the daylight environment that is routinely experienced despite the intervening dark periods (Dupuis & Ojeda et al. 2025). Here, we investigate how prior acclimation impacts algal fitness upon a change in daylight intensity and how quickly Chlamydomonas can reprogram its photoprotective strategy in a diurnal context. We performed a systems analysis of synchronized populations acclimated to diurnal LL when subjected to HL days and of populations acclimated to diurnal HL when subjected to LL days. In the latter case, diurnal photoacclimation decreased fitness during the first day at a new light intensity: HL-acclimated cells barely increased in size over the first LL period, and they failed to complete a cell cycle. However, although LL-acclimated cells showed severe photodamage after 6 hours of HL, they recovered chloroplast form and function later that afternoon and successfully divided at nightfall. These cells rapidly altered their thylakoid membrane ultrastructure, increased their photoprotective quenching capacity, and decreased their inventory of photosystem and antenna proteins by the end of the first HL day. Transcriptomic and proteomic analyses revealed rapid induction of thousands of genes, including those encoding proteases, chaperones, and other proteins involved in the chloroplast unfolded protein response. These results show that the alga is highly flexible and competent to rapidly acclimate to changes in diurnal light intensity.

All living cells depend on the dynamic balance between their energy supply and demand to survive and thrive in dynamic environments. In extreme cases, like for photosynthetic organisms, their energy source, light, can fluctuate dramatically in intensity over timescales of seconds to hours. While various photosynthetic electron flows (EF) are crucial for maintaining bioenergetic homeostasis, how EFs are modulated to respond to dynamic energy intake remains unclear. Here, we show in the model green alga Chlamydomonas reinhardtii that each EF is best suited to a specific domain of energetic fluctuation periodicity for which it can support the cells energetic needs, which we term bandwidth. By systematically exposing cells to a range of light periodicities, we show that while cyclic EF has a large bandwidth, pseudo-cyclic EF (PCEF) can only sustain the cells energetic needs for fast light fluctuations, and that the interplay between the chloroplast and the mitochondria (CMEF) has a limited bandwidth. We further show that the bandwidths of PCEF and CMEF, specialized for dynamic lights, are related to their capacity to generate ATP and protect the photosynthetic apparatus. Finally, we show that in wild-type cells, the activity level of PCEF and CMEF matches their bandwidth, and we propose that cells tune the relative activity of each AEF depending on the light fluctuation frequency. Our work opens an avenue of research to characterize the molecular mechanisms required for phototrophic growth in complex and dynamic energetic landscapes. It further provides a generalizable framework for understanding the physiological importance of molecular mechanisms in a dynamic environment.

Bacteria living on surfaces are often confined to droplets. When these droplets evaporate, the motion of the liquid-air interface and the associated internal capillary flow confine the bacteria. Here we study how E. coli bacteria interact with this capillary confinement and agglomerate at the droplets contact line. We identify three different types of bacterial pattern formation that depend on the bacterial activity and the environmental conditions imposed by the evaporating droplet. When the evaporation is fast, the bacteria are slow or the suspension is dilute, a uniform contact-line deposit forms. However, when the capillary confinement concentrates the bacteria at the contact line beyond a critical number density, localized collective motion spontaneously emerges. In that case, the bacteria induce a local stirring of the liquid that allows them to self-organize into periodic patterns and enables them to collectively escape from the contact line. At very high number densities, these periodic patterns get destabilized by bacterial turbulence in the bulk of the droplet resulting in the formation of mobile bacterial plumes at the contact line. Our results show how the subtle interplay between the bacteria and the capillary flow inside the droplet that surrounds them governs their dispersal. Significance StatementAn evaporating sessile droplet is a common natural habitat to bacteria. Bacteria that live inside the droplet are exposed to a confinement caused by the moving liquidair interface, and an evaporation-driven capillary flow that agglomerates them at the contact line. Here we show how bacteria interact with this confining flow. We identify three vastly different types of bacterial self-organization that depend on the bacterial activity and the environmental conditions imposed by the droplet. Our work is a first step towards understanding how the interplay between motile bacteria and the interfacial flows that exist in evaporating droplets affects their deposition onto surfaces, which is key to their future survival.

The reaction of CO2 with H2O to form HCO3- and H+ is one of the most important chemical equilibria in cells. In mammalian sperm, a soluble adenylyl cyclase (sAC) serves as cellular HCO3- sensor that conveys the equilibrium state via cAMP synthesis to cAMP-signaling molecules. The function of sAC and cAMP in non-mammalian sperm is largely unknown. Here, we identify sAC orthologs in sea urchin and salmon sperm that, surprisingly, are activated by alkaline pH rather than HCO3-. Two amino-acid residues required for HCO3- binding of mammalian sAC are lacking in pH-regulated sAC. Orthologs identified in ten other phyla are also lacking either one of these key residues, suggesting that pH control is widespread among non-mammalian metazoan. The pH-sensitive sAC controls several functions of sperm from external fertilizers. Upon spawning, alkalization triggers cAMP synthesis and, thereby, activates motility of quiescent sperm. Egg-derived chemoattractants also alkalize sperm and elevate cAMP, which then-modulates pacemaker HCN channels to trigger a chemotactic Ca2+ response. Finally, the sAC and the voltage- and cAMP-activated Na+/H+ exchanger sNHE mutually control each other. A picture of evolutionary significance is emerging: motility and sensory signaling of sperm from both internal and external fertilizers rely on cAMP, yet, their sAC is regulated by HCO3- or pHi, respectively. Acidification of aquatic habitats due to climate change may adversely affect pH-sensing by sAC and thereby sexual reproduction in the sea. Statement of significanceAdenylyl cyclases synthesize cAMP, a prominent cellular messenger. A bicarbonate-sensitive AC family member, soluble AC (sAC), is tied to the chemical equilibrium: H2O + CO2 {leftrightarrow} HCO3- (bicarbonate) + H+. The sAC is required for fertilization: Mammals lacking sAC are infertile and sperm immotile. We now identify a new sAC form in sperm of non-mammalian animals that reproduce in the sea. This novel sAC is activated at alkaline pH rather than bicarbonate. It controls sperm motility and chemotaxis. The switch from HCO3- to pH rests on substitution of two amino-acids, which represents an adaptation to aquatic environments low in bicarbonate. Acidification of aquatic habitats due to climate change may adversely affect sAC activity and, thereby, fertilization.

Regular light-dark cycles greatly affect organisms, and events like eclipses induce distinctive physiological and behavioural shifts. While well-documented in animals, plant behaviour during eclipses remains largely unexplored. Here we monitored multiple spruce trees to assess their individual and collective bioelectrical responses to a solar eclipse. Trees anticipated the eclipse, synchronising their individual bioelectrical behaviour hours in advance. Older trees displayed greater anticipatory behaviour with early time-asymmetry and entropy increases, characteristic of memory dynamics in open dissipative systems such as trees. These results reveal a complex relationship between individual and collective tree behaviour, shaped by individual age and physiology as well as collective history and memory. This highlights the adaptive significance of synchrony and collective memory in plants, offering new insights on coordinated behaviours in nature. One Sentence SummaryTrees anticipate and respond collectively to a solar eclipse.

Phosphate availability modulates plant immune function and regulates interactions with beneficial, phosphate-providing, microbes. Here, we describe the hijacking of plant phosphate sensing by a family of Nudix hydrolase effectors from pathogenic Magnaporthe oryzae and Colletotrichum fungi. Structural and enzymatic analyses of the Nudix effector family demonstrate that they selectively hydrolyze inositol pyrophosphates, a molecule used by plants to monitor phosphate status and regulate starvation responses. In M. oryzae, gene deletion and complementation experiments reveal that the enzymatic activity of a Nudix effector significantly contributes to pathogen virulence. Further, we show that this conserved effector family induces phosphate starvation signaling in plants. Our study elucidates a molecular mechanism, utilized by multiple phytopathogenic fungi, that manipulates the highly conserved plant phosphate sensing pathway to exacerbate disease. One-Sentence SummaryA family of conserved enzyme effectors from pathogenic fungi manipulate plant phosphate sensing to promote infection.

Cellular responses to abiotic stress involve multiple signals including secondary messengers such as reactive oxygen species (ROS) and Ca2+, phytohormones such as abscisic acid (ABA) and chloroplast-to-nucleus retrograde signals such as 3-phosphoadenosine 5-phosphate (PAP). Mechanism(s) by which these messengers, produced in different subcellular compartments, intersect for cell regulation remain enigmatic. Previously we showed that the chloroplast retrograde signal PAP, similar to ABA, induces an increase in ROS levels in guard cells (Pornsiriwong et al, 2017). Here we demonstrate a mechanistic link enabling ABA and PAP to coordinate both chloroplast and plasma membrane ROS production. In whole leaves, PAP alters various ROS-related processes including plasmodesmal permeability as well as responses to ozone and the bacterial elicitor flg22, but mainly initiates processes that quench ROS during oxidative stress. Conversely, we show in guard cells, both PAP and ABA induce an increase in ROS levels in both chloroplasts via photosynthetic electron transport, and the apoplast via the RESPIRATORY BURST OXIDASE HOMOLOG (RBOH). Both subcellular ROS sources were necessary for ABA- and PAP-mediated stomatal closure. However, PAP signaling diverges from ABA by activating RBOHD, instead of RBOHF, for apoplastic ROS mediated stomatal closure. We identified three calcium-dependent protein kinases (CPKs) as the post-translational activators of RBOHD-mediated ROS production. CPK13, CPK32, and CPK34 were transcriptionally induced by PAP and concurrently activate RBOHD and the slow anion channel SLAC1 by phosphorylating two Serine (S) residues, including S120 which is also targeted by the core ABA signaling kinase OPEN STOMATA 1 (OST1). Consequently, overexpression of the PAP-induced CPKs rescues stomatal closure in ost1. Our data identify stomatal chloroplasts, to be nodes in the multifaceted cellular stress response networks as they are both sources and mediators of ROS and retrograde signals such PAP. Thus, chloroplasts are not just mediators of photosynthesis in response to, for example, excess light, but can serve as critical nodes in the multifaceted cellular stress response networks in specialized cells via retrograde signals, providing support to the concept of sensory plastids. Significance StatementThe chloroplast is an environmental sensor for stresses such as excess light and drought via the activation of photosynthetic-mediated retrograde signals. However, how does it function in specialized cells for which carbon fixation is secondary? Here we show the chloroplast is an important node to coordinate multiple plant signaling pathways in response to stresses such as drought. The chloroplast retrograde signal 3-phosphoadenosine 5-phosphate (PAP) plays multiple roles in reactive oxygen species (ROS) signaling and homeostasis. While PAP suppresses ROS in photosynthetic tissue, PAP instead induces guard cell ROS in chloroplasts and extracellular space to induce stomatal closure. We decipher how PAP-induced proteins activate both extracellular ROS production and anion channels for stomatal closure, thus providing a mechanism by which chloroplasts provide a strategic complement to canonical hormonal pathways in regulating plant physiological responses in specialized cells.